Magnetization process for a quasi-one-dimensional S = 1 antiferromagnet 
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We investigate the magnetization process for a quasi-one-dimensional S = 1 antiferromagnet 
with bond alternation. By combining the density matrix renormalization group method with the 
interchain mean-field theory, we discuss how the interchain coupling affects the magnetization curve. 
It is found that the width of the magnetization plateau is considerably reduced upon introducing 
the interchain coupling. We obtain the phase diagram in a magnetic field. The effect of single-ion 
anisotropy is also addressed. 
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I. INTRODUCTION 

Low-dimensional quantum spin systems with the spin 
gap have attracted much interest recently. In particular, 
the Haldane system with integer spin@ has been stud- 
ied extensively. Q The origin of the Haldane spin-gap 
phase is now understood well in terms of the valence- 
bond-solid (VBS) picture. [g| Theoretical investigations 
have been extended@to more realistic systems, including 
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the single-ion 



e.g. the bond alternation 
the interchain coupling, 
anisotnw,0@,0 ) g,|l i 

In particular, the interchain coupling, which may trig- 
ger the quantum phase transition to a magnetically 
ordered phase, is indispensable to explain the ground 
state of some 5 = 1 compounds such as CsNiCl3,]T^] 
NENP,g§ PbNi 2 V 2 08,lPl SrNi 2 V 2 Q 8 .pl Theoreti- 
cally, Sakai and Takahashi]l3|] studied the effects of the 
interchain coupling on the Haldane system by combin- 
ing the exact diagonalization method with the interchain 
mean-field theory, and thus obtained the phase diagram 
for the quasi-one-dimensional (Q1D) system. More quan- 
titative treatments were later given by the series expan- 
sion method J8fl and also by the quantum Monte Carlo 
simulations. |J|7| 

Some striking phenomena for the quantum phase tran- 
sition found under an applied magnetic field have further 
stimulated the study on such Q1D spin systems. Most 
interesting may be the plateau formation in the magneti- 
zation curve at some fractional value of the full moment. 
Also, a long-range order induced by a magnetic field gives 
another interesting paradigm of the quantum phase tran- 
sition, which has been observed in some materials such as 
NDMAZ, @ NDMAP, Q NTEAP |22j and NTENP. 

In this paper, we study magnetic properties of a 
Q1D 5=1 Heisenberg antiferromagnet with bond al- 
ternation by combining the density matrix renormal- 
ization group (DMRG) with the interchain mean-field 
theory. jl3[ [l4|, gTj We clarify how the interchain coupling 
affects the magnetization curve with particular emphasis 
on the properties around the plateau. The effect of single- 
ion anisotropy is also discussed, which sometimes plays 



an essential role to understand low-energy properties of 
Q1D 5 = 1 compounds. 

This paper is organized as follows. In Sec. II, we in- 
troduce the model Hamiltonian for a Q1D 5=1 sys- 
tem, and briefly outline the numerical procedure. In 
Sec. Ill, we present the results, and discuss magnetic 
properties of the Q1D system. Furthermore the effect of 
single-ion anisotropy is discussed in Sec. IV. It is shown 
that the magnetization curve exhibits distinct behaviors 
depending on the relative direction between single-ion 
anisotropy and an applied magnetic field. Brief summary 
is given in Sec. V. 



II. 



MODEL AND NUMERICAL METHOD 



In order to study a Q1D 5 = 1 antiferromagnet in 
an applied magnetic field H, we consider the following 
Hamiltonian, 



n = j^2{i-(-i) i s}s jti -s j , i+1 
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where Sj i is the 5 



1 spin operator at the i-th site in 
the j-th chain. The summation with < j, j' > is taken 
over all the pairs of the nearest-neighbor chains. Here <5 
is the bond alternation parameter and J(J') the intra- 
chain (inter-chain) antiferromagnetic coupling. We set 
J = 1 and g/iB = L for convenience. 

We discuss magnetic properties of the above Q1D spin 
model by combining DMR G p5[ with the interchain 
mean-field theory. |l4] ^7$Tc 
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To deal with the effects 
of the interchain coupling in the presence of a mag- 
netic field, we introduce two kinds of mean fields as 
< 5f >^ (— l) l m s and < 5f >~ m u , where m s (m u ) 
is the staggered (uniform) moment induced by the inter- 
chain coupling and the applied magnetic field. We thus 
end up with the ID Hamiltonian, 



H 



MF 



2 



(-lfSf. (2) 



Here h s and h u are the effective internal fields, which 
are defined as h s = z^J'm s and h u = z±J'm u respec- 
tively, where zj_ is the number of the adjacent chains. 
The staggered (uniform) magnetization per site m s (m u ) 
is written down as 

m a = 4> s (h s ,h u ,S,H) , (3) 
m u = (j) u (h s ,h u ,S,H) , (4) 

where the functions 4> s and 4> u can be determined from 
the magnetization for the effective spin chain with both 
of the uniform and staggered fields. Since the mean-field 
Hamiltonian (^) is given as a function of H — h u , it is 
convenient to rewrite the staggered magnetization as 



m s = <j> s (h s ,8, H), 



(5) 



where the effective field is defined as H = H — h u . By 
solving the equation with the self-consistent condition 
h s = z±J'm s , we can discuss the quantum phase transi- 
tion from the spin-gap phase to a magnetically ordered 
phase. Namely, if this self-consistent equation has a fi- 
nite m s , a long-range order is induced by the interchain 
coupling. Otherwise, the staggered field induced by the 
interchain coupling becomes irrelevant, and thereby the 
spin-gap phase still persists against a magnetically or- 
dered phase. The critical value J' c for the phase transition 
is given in terms of the zero-field staggered susceptibility 
as, 



1 



z±J' 



dh s 



(6) 



This completes the formulation of the interchain mean- 
field treatment. 

We conclude this section by giving an explicit form of 
the staggered magnetization (f> s for the single chain case, 
which will play a key role in the following analysis. We 
calculate 4> s by means of the infinite DMRG method, 
which enables us to treat ID spin systems even when the 
total S z is not a conserved quantity. In Fig. |l|, we sum- 
marize characteristic profiles of the staggered magnetiza- 
tion (f> s with H being fixed. When H = 0, the system has 
the triplet-excitation gap over the dimer-singlet ground 
state stabilized by the bond alternation (6 = 0.5). In 
fact, the magnetization curve in staggered fields has a 
finite derivative at the origin, as seen from the solid line 
in Fig. implying that the dimerized ground state is 
stable up to a certain critical value. The introduction of 
an external field decreases the spin gap, and eventually 
triggers the quantum phase transition to a gapless phase 
(so-called Tomonaga-Luttinger phase) at H = 1.5, which 
is accompanied by the divergence of the spin susceptibil- 
ity shown as the dashed line in Fig. [j]. It is seen that the 
staggered susceptibility takes a finite value again in the 
large field region {H — 2.5), implying that the spin gap 
phase with the magnetization plateau is stabilized, for 
which a finite interchain coupling is necessary to destroy 
the spin-gap state. 
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FIG. 1: Staggered magnetization <f> 3 as a function of the stag- 
gered field h s for an isolated spin chain with bond alternation 
S = 0.5. The solid, broken and dot-dashed lines are the mag- 
netization curves for H = 0.0, 1.5 and 2.5, respectively. 



III. EFFECTS OF THE INTERCHAIN 
COUPLING 



We now discuss how the interchain coupling affects 
the magnetization process in the Q1D system, follow- 
ing the procedure outlined above. In the following we 
will explicitly use the bare external field H instead of H 
(H = H + h u ). 



In Fig. ^ we show the magnetization curve calculated 
for the Q1D system, which is compared with that for the 
pure ID case (J' = 0). Note that the ground-state for the 
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FIG. 2: Magnetization curves for 5 — 0.5 as a function 
of the applied field H. The solid and broken lines are the 
results calculated for an isolated chain (z±J' = 0) and weakly 
coupled chains (z±J' = 0.04), respectively. 

ID case is either in the singlet-dimer phase (S > 0.25) or 
in the Haldane phase (6 < 0.25). We show the magnetiza- 
tion for the dimer-singlet phase (5 = 0.5) as an example, 
since that for the Haldane phase exhibits a similar behav- 
ior. Let us begin with the magnetization for the ID case. 
Since the system is in the singlet-dimer phase with the 
spin excitation gap, the uniform magnetization does not 
appear up to the critical field H c q. Beyond the critical 
field, the magnetization begins to increase and the sys- 
tem enters the Tomonaga-Luttinger liquid phase, which 
is immediately driven to a magnetically ordered phase 
upon introducing the interchain coupling. Further an in- 
crease in the field drives the system to another plateau 
state at m u = 1/2. This spin-gap phase with m u = 1/2 
is expected to be stable against the introduction of the 
interchain coupling. However, it should be noticed that 
the spin-gap vanishes at both ends of the magnetization 
plateau, so that the width of the plateau should be con- 
siderably reduced by the interchain coupling. We can see 
that this is indeed the case for the results in Fig.|[ 

To see the above behavior in more detail, we show 
the phase diagram for the QID spin system with the in- 
terchain coupling z±J' = and 0.04 in Fig. |[ Note 
that our phase diagram for the isolated chain (z± J' = 0) 
agrees fairly well with that already obtained by the ex- 
act diagonalization method. |ll[] Let us consider how 
the interchain coupling affects the plateau region with 
m u = 1/2, which is stabilized by the bond alternation. 
As seen from Fig. ||, the region of this phase is extended 
with increasing S. On the other hand, when the inter- 
chain coupling is introduced, the quantum fluctuation 
stabilizing the plateau is suppressed, and thus the region 
for the m u = 1/2 plateau is reduced. 

We note here that the phase diagram obtained here 
is qualitatively valid except for an extreme case, 8 ~ 1, 
where the system is approximately described by an as- 
sembly of isolated dimers for J' = 0. When a finite inter- 
chain coupling J' is introduced for such almost isolated 
dimers, the system first forms the ladders rather than 



FIG. 3: Phase diagram for the QID spin system. The solid 
and broken lines represent the phase diagram for an isolated 
chain (z±J' = 0) and weakly coupled chains (z±J' = 0.04), 
respectively. The uniform magnetization takes the quantized 
value at m u = in the region from to H c \ {H' cl ), m u = 1/2 
in the region from H C 2 {H' c2 ) to H C 3 (H' c3 ), and m u = 1 
above H„ {H' g ). In other regions, the magnetization varies 
continuously. 



the chains for two-dimension. In this case, our approach 
based on the interchain mean-field theory may not be 
sensible. To treat the critical phenomena in this region, 
we need more precise treatment beyond mean-field the- 
ory, which is now under consideration. 

Before closing this section, a brief comment is in order 
for the magnetization process in the QID Haldane system 
with a ferromagnetic interchain coupling, J' < 0. By 
rewriting the effective Hamiltonian, the critical points of 
the system with the ferromagnetic interchain couplings 
are given as, 

H' cQ {J' = -\J'\) = H' cQ {J' = \J'\) (7) 

H' c xA J> = "I J 'D = H 'ciA J ' = l J 'D - ^1 J '\ ( 8 ) 
H' s (J' = -\J'\) = H' s (J' = \J'\)-2z x \J'\. (9) 

As easily seen from these formulae, the width of the 
plateau H' c2 — H' cl does not depend on the sign of J' 
in the framework of our mean- field treatment. 



IV. EFFECTS OF SINGLE-ION ANISOTROPY 

In this section, we discuss the effect of single-ion 
anisotropy with the Hamiltonian, Hd — D Ylii^i) 2 - We 
consider two cases of anisotropy with respect to the direc- 
tion of an applied magnetic field, i.e. D _L H (a = y) and 
D || H (a = z). By choosing the x-axis as the axis for 
the magnetic order, we discuss the magnetic properties 
of the QID system. 

To begin with, let us consider a simple system in the 
isolated-dimer limit with (S,J f ) — (1,0), which provides 
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FIG. 4: Magnetization curves for the isolated dimer system 
{6 = 1). For D || H(D = 0.5), D _L H(D = 0.5) and D = 0. 
The inset shows the behavior around the plateau m u — 1/2. 



us with some insight into the role of single-ion anisotropy. 
The magnetization curves for D = 0.5 are shown as the 
broken line (D \\ H) and the solid line (D J_ H) in Fig^. 
Note that the magnetization shows jumps in a magnetic 
field, reflecting the isolated dimer-singlet ground state 
When D \\ H, the plateaus with m u — 0, 1/2 and 1 
appear in the magnetization process clearly. On the other 
hand, in the case of D _L H, the plateau is not exactly 
flat but is smeared by the Van Vleck contribution, as seen 
from the enlarged picture in the inset. This is because the 
total S z is not a conserved quantity for the case of D _L H 
even for h s = 0, which is in contrast to the case of D \\ H. 
We also find that the width of the m u — 1/2 plateau is 
increased (decreased) upon introducing anisotropy, D \\ 
H (D _L H). Therefore, the observation of the plateau 
structure may be somewhat difficult for D _L H in generic 
Q1D spin systems, because of shrinking of the width as 
well as smearing due to the Van Vleck term. 

Keeping the above properties in mind, we now consider 
the magnetization process for the pure ID as well as the 
Q1D systems. The numerical results are summarized in 
Fig. U Shown in Fig. |(a) are the magnetization curves 
for D || H(D = 0.5) and S = 0.5, Since the presence 
of anisotropy, D || H, has a tendency to stabilize the 
m u — 1/2 plateau, as mentioned above, thus increasing 
the plateau width compared with the isotropic case D = 
(see Fig. |^). It is also seen that the plateau state is 
rather stable against the introduction of the interchain 
coupling. 

In Fig. H (b), we show the magnetization curves for 
the case of D J_ H with bond alternation 8 = 0.7. In 
contrast to the above D \\ H case, the plateau is indeed 
obscured both for the ID (solid line) and Q1D systems 
(dashed line), which is due to the Van Vleck contribu- 
tion. Therefore, in order to obtain the phase transition 
point correctly, it is necessary to consider another proper 
quantity characterizing the transition clearly. For this 
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FIG. 5: Magnetization curves for the Q1D system: (a) D || 
H(D = 0.5) and 5 = 0.5, and (b) D _L H(D = 0.5) and 
S = 0.7. The solid and broken lines show the results for 
z±J' — and 0.04, respectively. 



purpose, we evaluate the staggered magnetization in the 
x-axis. The calculated results for a single-chain system 
are plotted (broken line) together with the uniform mag- 
netization (solid line) in Fig. |[ We note that a negligibly 
small staggered field h s — 1.0 x 10~ 5 is added to stabilize 
the spontaneous staggered magnetization. It is seen that 
although the plateaus at m u = 0, 1/2 and 1 are vague, 
the phase transition point can be determined clearly at 
which the staggered magnetization vanishes. This is used 
to determine the phase diagram for the case of D _L H. 
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FIG. 6: Uniform and staggered magnetizations for an iso- 
lated spin chain when D _L H (D = 0.5). 

Finally, we show the 5-H phase diagram for D \\ H and 
D _L H in Fig. ^ (a) and (b). In these figures, the solid 
and broken lines give the phase diagrams for an isolated 
ID chain and weakly coupled chains, respectively. We 
note that when J 1 = H = 0, the quantum phase tran- 
sition between the Haldane phase and the singlet-dimcr 
phase occurs at the critical point 5 C ~ 0.23 for D = 0.5. 
0, [[I In the case of D \\ H, the effects of the inter- 
chain couplings are analogous to those for the D = 
case (see Fig. ||). On the other hand, when the system 
has anisotropy in the y-direction, i.e. D _L H, the phase 
diagram exhibits a quite different feature in contrast to 
the D || H case. In particular, the m u = 1/2 plateau 
region is reduced considerably both for the pure ID as 
well as the Q1D systems. 

V. SUMMARY 

We have studied the magnetization process for a Q1D 
5 = 1 antiferromagnet with bond alternation. By com- 
bining the density matrix renormalization group method 
with the interchain mean-field treatment, the magnetic 
properties have been discussed. It has been shown that 
the introduction of the interchain coupling enhances 
the antiferromagnetic correlation, thereby reducing the 
width of the plateau in the magnetization curve. 

We have also discussed the effects of single-ion 
anisotropy, which sometimes plays an important role 
for understanding experimental findings in S = 1 spin 
systems. It has been clarified that the magnetization 
curves exhibit distinct behaviors depending on the rela- 



tive direction between single-ion anisotropy and an ap- 
plied field. In particular, we have found that the magne- 
tization plateau in the case of D _L H is rather unstable 
against the interchain coupling, and also suffers from the 
smearing effect due to the Van Vleck contribution, mak- 
ing the experimental observation somewhat difficult. 
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ci 1 — 1 1 1 1 1 1 1 
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FIG. 7: Phase diagrams for the Q1D quantum spin chain 
with single-ion anisotropy. The solid and broken lines are the 
phase boundaries for an isolated chain (z±J' = 0) and weakly 
coupled chains (z±J = 0.04), respectively. The meanings of 
the critical fields and the corresponding phases are as in Fig. 
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